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DECUS Program Library Write-up 


DECUS No. 8-293 


Coordinate Program ' 

This program consists of a main program and a large number of 
subroutines which are designed to facilitate developing the coordinates 
of the atoms of a given molecule. Loading information is allowed by a 
description of the available options, and of the algorithms which are 
used. The program occupies essentially all of core except for that 
used by version C of the floating point interpreter. Some locations 
normally used by the interpreter are also used, so that it is necessary 
to load the interpreter before loading the program. The binary tape which 
is supplied contains both the interpreter and the program. 

Loading 

The binary tape is loaded in the usual fashion via the binary loader. 
The starting address is 200. In case of an error, recovery can be made 
by restarting at 4746 or 1026. If this does not work, reload the program 
(this does not erase data) and restart at one of these locations. 

Operation and Options 

On starting, the computer will type 

ATOM 1 IS AT 0, 0, 0 
ATOM 2, R(l, 2) = 

The 1-2 bond length is then typed followed by a space. The computer now 
types 

ATOM 3, R(2, 3) = 

❖ The development of this program was assisted by a National Science 


Foundation Grant. 










The bond length is typed followed by a space. The response is 
ANGLE (1, 2, 3)=. The bond angle is typed. Atoms 1, 2 and 3 are 
placed in the x, y plane, with atom 2 lying on the x axis. Further atoms 
require information on bond length, bond angle and torsional angles. 

Thus, the computer types 

NOW TYPE R(K, L), ANGLE(J, K, L) and TORS (I, J, K, L) WHERE L IS 
THE NEW ATOM 
I J K L R(K, L) ANGLE TORS 

Integer values for the indices are now typed followed by decimal values 

for the bond length, angle and torsional angle. Each value is followed by 

a space. As an example, for cyclohexane with tetrahedral angles, the 

initial values of R(l, 2)=1. 54, R(2, 3)=1. 54 and ANGLE(1, 2, 3)=109. 5 

would be followed by 

01 02 03 04 1. 54 109. 5 60. 0 

02 03 04 05 1. 54 109. 5 -60.0 

03 04 05 06 1. 54 109. 5 60. 0 

The coordinates of the carbons will now have been developed. An exit 

from this mode of operation is effected by giving the first index a value 

of -1. The computer then types > and waits for the control letter followed 

by a space. The set of control letters and their effect is given in Table I. 

As a check on the coordinates defined by the values given above (note 

that one torsional angle must be negative in order to close the ring), we 
b b.b 

may type R 1 6 . (The superscript b indicates a space). The computer 
then answers with 1. 54 indicating the correct distance. One bond has not 
been defined by the data given above - that between atoms 1 and 6. This 
information is needed if the coordinates are to be plotted later. To add 
this type E, followed by a space. The computer responds with 
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ATOMS 

00010002 

00020003 

00030004 

00040005 

00050006 

ADD 

b b 

Now, type 1 6 to insert the additional bond. Terminate mode by 

typing -1 followed by a blank. The computer responds with DELETE. 

Since no pair of atoms is to be deleted, respond with -1 followed by a 

blank. The computer now types > and waits for another command. 

The secondary hydrogens may now be added. The sequence would be 

S 6 1 2 1. 1 106. 0 

S 1 2 3 1. 1 106. 0 

S 2 3 4 1. 1 106. 0 

S 3 4 5 1. 1 106. 0 

S 4 5 6 1. 1 106. 0 

S 5 6 1 1.1 106.0 

Here, S indicates that two hydrogens are to be added to the second atom of 

the list that follows (cf Table I). The three numbers give the indices of the 

atoms forming the C-C-C bond to which the hydrogens are to be added. The 

next number gives the bond length and is followed by the H-C-H bond angle. 

The secondary hydrogens will now have been added using the specified 

bond length and angle. The set of coordinates which have been developed 

b b 

may be outputted by typing P 0 (zero). If the coordinates are to be 
punched on tape for later tape to card conversion, type The 

computer halts, and the punch may be turned on. Pressing CONTINUE 
will cause the coordinates to be printed and punched in decimal form, and 
the computer again halts. If the tape is to be loaded back into the computer 
at a later time, shut off punch, put teletype on local, turn on punch and type 
00 . In either case, the program is continued by pressing CONTINUE. 
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The other options should be fairly clear from the descriptions in 

Table I. The one designated by L is particularly useful. If one wishes to 

replace a hydrogen with a methyl group, for example, one may translate 

so that the coordinates of the carbon to which it is attached is 0, 0, 0. Then 

L followed by the coordinates of the hydrogen will rotate the entire set of 

coordinates so that the hydrogen will lie on the z axis. Using the change 

(C) option, the hydrogen may be replaced by a carbon with the appropriate 

C-C bond length. The coordinates of methyl hydrogens may be added by 

translating to make the coordinates of the new carbon 0, 0, 0 and adding 

the following using the fix (F) option: 

01 0. 0 -1. 0343 0. 374 
02 0. 8957 0. 5172 0. 374 
03 -0. 8957 0. 5172 0. 374 
00 

Note that a zero followed by a number and a blank must precede each set 
of floating point numbers. This is necessary so that the computer may 
ignore leader, line feed and carriage return characters when a tape is 
read in using the F option. The computer initially checks each characte 
to see if it is a zero. When it finds the first zero, it reads the following 
number. If the number is greater than zero, the x, y and z coordinates 
are read using the floating point interpreter. If the number is zero, an 
exit from the subroutine is called and another option character is requested. 
If a particular orientation of the methyl group with respect to the rest of the 
moleciile is desired, the molecule may be suitably rotated about the z axis 
(U b 3 b v where v is the number of degress) before adding the hydrogen 
coordinates. 


4 








The program will accept up to 42 atoms. A larger number will result 
m destroying data and part of the program. No check is provided to guard 
against this possibility. Particular attention should be paid to the bond 
register. If a large number of changes in coordinates are made, it is 
possible for this register to overflow without exceeding 42 atoms. It 
should therefore be checked occasionally - and it will have to be correct 
if the WRITE program is to operate correctly. In deleting coordinates 
using the delete (D) option, it is important to note that an attempted 
deletion of non-existent data (i. e. n t > no. of atoms) results in deleting 
the coordinates of the last atom in the list. 

It is possible to have difficulty entering hydrogens using the S 
option for certain arrangements of atoms. This is not a frequent 
occurrence and may be corrected by deleting the incorrect atoms, rotating 
the coordinates and using S again. The bond register also will have to be 
corrected. 

The WRITE (on CRT) option is described in a separate writeup. 

A sample run using this program is appended . 
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CYCLOHEXANE 




On 


be-fv 

in i ti ai 


looxl before starting program 


CALCULATE ATOMIC COORDINATES 
ATOM 1 IS AT 0,0,0 
ATOM 2, R(l,2) = 1.54 
ATOM 3, R<2,3> = 1.54 ANGLE (1,2,3) 
NOW TYPE R(K,L>, 

THE NEW ATOM 


procjCa m 


ft yj i*eoJ\ 


=109.5 


ANGLE(J,K,L) AND TORSCI, J,K, L) WHERE L IS 


J 


i S+ 


three aTOmS 
are defined 


I 

J 

K 

L 

RCK,L> ANGLE TOR: 

01 

02 

03 

04 

1.54 

109.5 

60*0 

02 

03 

04 

05 

1.54 

109.5 

-60*0 

03 

04 

05 

06 

1.54 

109.5 

60»0 


il£/"£ y the Carbon skeleton i$ built up 


-1 


> R 1 6 +0.1542879E+01 tf^ks for bond l&wjfb (undefined bond IS checked) 

> A 2 1 6 +0.1094618E+03 A asVs for bond am<^le [undefined crnjle is ebeebed) 
— £ (XL ks for iishnj of bond register 


>E —s- 

ATOMS 

00010002 

00020003 

00030004 

00040005 

00050006 

ADD 1_6 

ADD M 
DELETE - 


1 


>s 

6 

1 

2 1.1 

106 

>s 

1 

2 

3 1.1 

106 

>s 

2 

3 

4 1.1 

106 

>s 

3 

4 

5 1.1 

106 

>s 

4 

5 

6 1 

i.l 

106 

>s 

5 

6 

1 1 

l.l 

106 


atoms cob/d ace bonded j ie. 1-Z Z-5 etc. 


add the I jlo bond 
-l derm mates mode 

-I -termmates mode 


) 


^ pu4"S In SecOnd-OCoj h^jdcooje 1m£ on middle afom O' 

liU ; cJltb I I & bond ItMOpW and /66° bond ounces 


Qsncjl 


> P 0 
ATOM 
0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 


J 


J 


+0.0000000E+00 
+0.1539999E+01 
+0.2054054E+01 
+0.1541440E+01 
+0.1441717E-02 
-0.5140514E+00 
-0.3815353E+00 
-0.3815352E+00 
+0.1922070E+01 
+0.1922070E+01 
+0.1719231E+01 
+0.3153556E+01 
+0.1924548E+01 
+0.1923544E+01 
-0.3805962E+00 
-0.3795888E+00 
-0.1613584E+01 
-0.1802634E+00 


+0•0000000E+00 
+0.0000000E+00 
+0.1451670E+01 
+0.2178525E+01 
+0.2180052E+01 
+0.7288928E+00 
+0.4892026E+00 
-0.1031332E+01 
-0.5406136E+00 
-0.5406136E+00 
+0.1975553E+01 
+0.1467640E+01 
+0•3209266E+01 
+0.1688022E+01 
+0.2720688E+01 
+0.2721399E+01 
+0.7146030E+00 
+0•2059813E+00 


+0.0000000E+00 
+0.0000000E+00 
+0.0000000E+00 
-0.1257180E+01 
-0.1258062E+01 
-0.1258945E+01 
+0.9083559E+00 
+0.2801060E-01 
+0.8784978E+00 
-0•8784978E+00 
+0.9074350E+00 
+0.2893292E-01 
-0.1285477E+01 
-0.2164595E+01 
-0•3795649E+00 
-0.2136560E+01 
-0.1287592E+01 
-0.2167322E+01 


huwber rjnces J u sf 
■fbocf coordinate / and 
-/ cjt/eS (x dt-ffeYOnf 
formed printout tollaJi 
a HALT (press Cout) 
do aMous puncL to be 
du/ned om 


> Another command is requested 
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Table I 


5jC 

Control Characters and Their Effect 


Character and List 


Effect 


R 

n l 

«2 


A 

n i 

n 2 

n 3 

T 

n l 

n 2 

n 3 n4 


S nj n 2 n 3 vj v 2 

M ni n 2 n 3 n 4 

G 

C n x V! v 2 v 3 

D n x 

I n x v 2 v 2 v 3 

F nj V! v 2 v 3 


The distance between atoms n x and n 2 
is typed. 

The angle formed by atoms nj, n 2 and 
n 3 is typed. 

The torsional angles formed by atoms 
n x , n 2 , n 3 and is typed. 

Secondary hydrogens are added to atom 
n 2 which is attached to n x and n 3 . The 
bond length is v* and the HCH bond angle 
is v 2 . 

A methine hydrogen is added to atom n x 
which is attached to atoms n 2 , n 3 and n 4 . 

The bond length is v*. 

Get additional coordinates. A new set 
of indices, bond length, bond angle and 
torsional angle is then typed. Exit from 
this mode is made by typing -1 as before. 
Change the coordinates of atom n x to v^, 
v 2 , v 3 . 

Delete the coordinates of atom n x . If n x 
is greater than the number of atoms, the 
coordinates of the last atom are deleted. 
Insert the coordinates vj, v 2 , v 3 before 
those of atom n 1# 

Fix the coordinates of another atom as 
v i> v 2 ’ v 3 * Additional coordinates may be 
added until a zero value of n x is reached. 
The computer then asks for a new in¬ 
struction. n x must be preceded by a zero 
in each case. Thus, the last entry is 00 
followed by a space. 

Examine bond registers. The attached 
atoms are listed in pairs. The computer 
then types ADD. Additional bonding pairs 
may now be entered. Adding is terminated 
by -1. The computer then types DELETE. 
The number in the list of the entry to be 
deleted may now be typed. This mode is 
terminated by -1. 
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p nl 

Print coordinates. If n x is zero, all 
coordinates are typed in floating point 
form. If nj is positive, the coordinates 
of atom nj are typed. If n x is negative, 
the computer halts, permitting the punch 
to be turned on. Pressing the CONTINUE 
key leads to typing of coordinate values 
in decimal form. The computer again 
halts and is restarted by pressing the 
CONTINUE key. If the tape is to be read 
back into the computer at a later time, 
terminate the tape with 00 (see text). 

w 

Write coordinates in CRT. This leads 
to the message ’’LOAD SECOND TAPE, 
START AT 4400. ” See writeup of the 
WRITE subroutine. 

B n| 

Invert the values of the x(n x = l), y(nj=2) 
or z(nj=3) array (i. e. multiply them by 
-1.0). 

Li Vi v 2 v 3 

Rotate the vector defined by vj, v 2 , v 3 so 
that it lies along the z axis. 

Q nj n 2 

Transpose arrays nj and n 2 (i. e. ex¬ 
change them). For x, n=l; y, n=2; 
z, n=3. 

U rij vj 

V n, vi v 2 v 3 

Rotate about axis nj by v* degrees. 

If n x is a positive integer, translate 
coordinates of atom n to 0, 0, 0- 
ignore Vj, v 2 , v 3 . If nj is zero or negative, 
translate coordinates to make v x , v 2 , 
v 3 =0, 0, 0. 


* 

The values designated as n are integers, those designated by v are decimal 
numbers, A space or carriage return must follow each entry. 

If the control character is entered incorrectly, type it a second time. The 
computer will respond with a 11 ? M , followed by > . Then a new character 

may be entered. 

If an integer is entered incorrectly, type a back arrow («-) before the 
terminating space. Then type the entire correct number. If the space is 
typed, stop computer and restart at 4746. 

If a decimal (floating point) number is entered incorrectly, type a ’’rubout” 
before the terminating space. Then type the entire correct number. 


8 






Method of Calculation 


The first atom is located at 0, 0, 0 and its bond to the second 
atom is considered to lie on the x axis. The third atom is placed in the 
xy plane with x 3 = x 2 - r 23 cos0 and y 3 = sin0. 

The remaining atoms are located as follows where R is the 
bond length between atoms k and 1, 0 is the bond angle formed by atoms 
j, k and 1, and 4> is the torsional angle formed by atoms i, j, k and 1. 

A new coordinate system is set up so that 

x^ ! = Rcos0 
y^'= Rsin0cosc|) 
z^*= Rsin0sinc(> 

Then the transformation matrix is obtained which will transform the 
new coordinate system into the original system. Thus: 


X 1 

= X. 


T • 

V 

+ 

T • 

y « 

+ T • 

z. 1 

1 

k 


XX 

1 


xy 


xz 

1 

y. 

= Yi 

+ 

T • 

x l* 

+ 

T • 

y • 

+ T • 

z ! 

1 

k 


yx 

1 


yy 

■^1 

yz 

1 

z. 

= z n 

+ 

T • 

x l ! 

+ 

T • 

y* 

+ T • 

z * 

1 

k 


zx 

1 


zy 

^1 

zz 

1 


The transformation matrix is obtained as follows: 


R ab ='^y\F‘+ <y j -y 1 7 T <y v T 


L* - <Y x k>/ R ab 


T =(y.-y)/R 
yx j k ' ab 


T ZX = <yV/ R ab 


SCALE = (x.-x, )T + (y.-y, )T + (z. -z )T 

l k xx w i 7 k yx l k zx 


R, = (x.-x, ) - SCALE- T 
1 l k xx 

R 2 = (y.-y, ) - SCALE - T 

i k yx 


r 3 = ( z i-z k ) - SCALE - T 


zx 
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RAJ = + R 2 2 '+ R 3 t ~ 


T 


xy 


Rj/ RAJ 


T = R J RAJ 

yy 


T = Rj/ RAJ 
zy 



T • T 


zx 


yy 


rp _ p • rp 

yz zx xy 


T • T 


XX 


zy 


T = T • T 


- T * T 

yy yx xy 


ZZ XX 


The bond lengths and bond angles are calculated using the usual 
relationships. The dihedral angle (<|>) is calculated as follows where atoms 
1,2, 3and 4 define the angle: 


c = [ (yi-y 2 )(z 2 - z 3)-( z i- z 2)(y2-y3)]* [ (y4-y3)(z2-z 3 )-(z 4 -z3)(y 2 -y 3 )] 

+ [ ( Z 1 “ z 2 ) (^2 -^3 ) “ (x t ”X 2 ) (z 2 - Zj ) ] * [ (z 4 -Z 3 )(x 2 -X 3 )-(x 4 -X 3 )(z 2 -Z 3 )] 

+ [ (x 1 -x 2 )(y 2 -y3)-(y 1 -y 2 )(x 2 -x 3 )]- [ (x 4 -X3)(y 2 -y 3 )-(y 4 -y3)(x 2 -X3)] 

D = {[ (yi-y2)( z 2- z 3)-( z l-z 2 )(y 2 -y3)] 2 + [ ( z 1 -Z 2 Hx 2 -X 3 )-(x 1 -X 2 )(z 2 -Z 3 )] 2 
* 1 

+ [ (x a -x 2 )(y 2 -y 3 )-(y 1 -y 2 )(x 2 -X3 )] 2 } 2 • {[ (y 4 -y 3 )(z 2 -Z3)-(z 4 -z3)(y 2 -y 3 )] 2 
+ [ (z 4 -Z3)(x 2 -x 3 )-(x 4 -x 3 )(z 2 -Z3)] 2 + [ (x 4 -X3)(y 2 -y3)-(y 4 -y 3 )(x 2 -X3)] 2 } 2 


<|> = cos” 1 (C/ D) 
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Coordinate 


0002-0037") 

0063-0177J 

0200-1023 

1025-1177 

1200-1265 

1266-1277 

1300-1327 

1330-1372 

1373-1457 

1460-1562 

1563-1633 

1634-1732 

1733-1747 

1750-1777 

2000-2060 

2061-2357 

2360-2574 

2600-2641 

2642-2674 

2675-2777") 

4400-4520J 

4521-4606 

4607-4722 

4723-4745 

4746-4750 

4315-4376] 

4143-4176J 

3654-3703 

3704-3777 

5000-5123) 

5252-5376J 

5544-5562 

4751-4756 

3000-3177 

3200-3377 

3400-3577 

3600-3652 

4000-4057 

4075-4142 

4200-4311 

4757-7577 


Program Storage Map 


Page 0 constants 


Initial operating program 


PRNT- 

ANGLES- 

RVAL- 

SDIST- 

TRPOSE- 

DELETE- 

ROT2- 

TRANSL- 

LOCATN- 

CHANGE- 

FIX- 

PRINTT- 

SECOND- 

TORSS- 

LIST- 

ARCOS- 

ROT1 - 

INSERT - 
TERT- 
WRITE- 
RECOVR- 

EX AMINE 

INVERT - 
PRNT2- 

DCOUT- 

SIGN- 

ERR- 


Subroutine control program 
Bond angle subroutine 
Bond distance subroutine 

Subroutine for calculating distance between points 
Subroutine for transposing coordinates 
Subroutine for deleting coordinates 
Subroutine for rotating about an axis 
Subroutine for translating coordinates 
Subroutine for obtaining addresses of elements 
in x, y and z arrays 
Subroutine for changing a coordinate 
Subroutine for fixing additional coordinates 
Subroutine for printing coordinates using E format 
Subroutine for adding secondary hydrogens 
Torsional angle subroutine 

Subroutine which prints coordinates - used with PRINTT 
ARC COSINE subroutine 

Subroutine for rotating a point onto the Z axis 

Subroutine for inserting a coordinate 
Subroutine for introducing tertiary hydrogens 
Dummy write subroutine (see text) 

Error recovery subroutine 

Subroutine for examining bond register 

Subroutine for changing sign along one coordinate 
Print selection subroutine 

Decimal output subroutine 

Subroutine used by ROT1 
Subroutine used by DC OUT 


X Array 
Y Array 
Z Array 
Bond register 

MESAGE- Typeout subroutine (DIGITAL-8-18) 

DECPRT- Decimal print subroutine (DIGITAL-8-22) 

SICONV- Decimal input subroutine (DIGITAL-8-28) 

Floating point interpreter »C" (with overlays from some of 
the above subroutines) 
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0002 

0003 

0004 

0005 

0006 

0007 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0030 

0031 

0032 

0033 

0034 

0035 


/PROGRAM FOR OBTAINING CARTESIAN COORDINATES 

CALL= 4400 

RETURN=5400 

F SQRT=2 

F CO S= 4 

FSI N= 3 

F ATAN= 5 

BAX=146 

BAY=151 

BAZ=154 

TZ 1 = 65 

TZ2=157 

TZ3=143 

INPUT=13 

OUTPUTS 14 


FLPT= 4407 
*2 



3000 LOCX, 
3200 LOCY, 
3400 LOCZ, 
3600 LOCR, 


3000 

3200 

3400 

3600 

0 

5600 


0000 NN, 
5600 


0000 LX, 
0000 LY, 
0000 LZ, 
0000 IX, 
0000 IY, 
0000 IZ, 
0000 JX, 
0000 JY, 
0000 JZ, 
0000 KX, 
0000 KY, 
0000 KZ, 
0000 IR, 


0 

0 

FLTG 0.0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 



0000 ZERO, 


0000 

0000 


0000 R, 
0000 


FLTG 0.0 


0000 

7773 CONST, 


FLTG 0.01745329 


2167 

6416 
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0036 

4000 

MESAGE, 

4000 


0037 

0000 

N, 

0 




*63 



0063 

0000 

INDEX, 

0 


0064 

0003 

THREE, 

3 


0065 

0000 

ANGLE, 

FLTG 

0 • 0 

0066 

0000 




0067 

0000 




0070 

0000 

TORS, 

FLTG 

0.0 

007 1 

0000 




0072 

0000 




0073 

0000 

TEMP, 

FLTG 

0.0 

0074 

0000 




0075 

0000 




0076 

0000 

II, 

0 


0077 

0000 

J, 

0 


0100 

0000 

K, 

0 


0101 

0000 

L, 

0 


0102 

4200 

SICONV, 

4200 


0103 

407 5 

DECPRT, 

4075 


0104 

0000 

FLAG, 

0 


0105 

0000 

TX1, 

FLTG 

0.0 

0106 

0000 




0 107 

0000 




0110 

0000 

TX2, 

FLTG 

0.0 

0111 

0000 




0 112 

0000 




0113 

0000 

TX3, 

FLTG 

0.0 

0114 

0000 




0115 

0000 




0116 

0000 

TY1, 

FLTG 

0.0 

0117 

0000 




0120 

0000 




0121 

0000 

TY2, 

FLTG 

0.0 

0 122 

0000 




0 123 

0000 




0124 

0000 

TY 3, 

FLTG 

0.0 

0125 

0000 




0126 

0000 




0127 

0000 

SS, 

FLTG 

0.0 

0130 

0000 




0131 

0000 




0132 

0000 

CS, 

FLTG 

0.0 

0133 

0000 




0134 

0000 




0135 

0000 

C2, 

FLTG 

0.0 

0136 

0000 




0137 

0000 




0140 

0000 

RAJ, 

FLTG 

0.0 

0141 

0000 




0142 

0000 




0143 

0000 

RAB, 

FLTG 

0.0 

0144 

0000 




0145 

0000 








0146 

0000 

JAX * 

FLTG 

0*0 

0147 

0000 




0150 

0000 




0151 

0000 

JAY * 

FLTG 

0.0 

0152 

0000 




0153 

0000 




0154 

0000 

JAZ* 

FLTG 

0.0 

0155 

0000 




0156 

0000 




0157 

0000 

SCALE* 

FLTG 

0.0 

0 160 

0000 




0161 

0000 




0162 

0000 

SUM* 

FLTG 

0.0 

0163 

0000 




0164 

0000 




0165 

0000 

CAX* 

FLTG 

0.0 

0166 

0000 




0167 

0000 




0170 

0000 

CAY* 

FLTG 

0.0 

0171 

0000 




0172 

0000 




0173 

0000 

CAZ* 

FLTG 

0.0 

0174 

0000 




0175 

0000 




0176 

0513 

CONT* 

COORD 

0177 

1025 

PRINT* 

PRNT 




*3600 



3600 

0102 


0102 


3601 

0203 


0203 







